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ABSTRACT 

This  paper  reports  new  Very  Large  Array  (VLA)  and  Very  Long  Baseline  Array  (VLB  A)  observations  of  the 
extragalactic  source  B1 849+005  (Galactic  coordinates  i  =  33?44,  b  =  +0?21)  at  frequencies  between  0.33 
and  15  GHz  and  the  reanalysis  of  archival  VLA  observations  at  0.33,  1 .5,  and  4.9  GHz.  The  structure  of  this  source 
is  complex,  confirming  previous  suggestions,  but  interstellar  scattering  dominates  the  structure  of  the  central 
component  at  least  to  15  GHz.  An  analysis  of  the  phase  structure  functions  of  the  interferometric  visibilities  shows 
the  density  fluctuations  along  this  line  of  sight  to  be  anisotropic  (axial  ratio  =  1.3)  with  a  frequency-independent 
position  angle  and  having  an  inner  scale  of  roughly  a  few  hundred  kilometers.  The  anisotropies  occur  on  length 
scales  of  order  1015  cm  (D/5  kpc),  which  within  the  context  of  certain  magnetohydrodynamic  turbulence  theories 
indicates  the  length  scale  on  which  the  kinetic  and  magnetic  energy  densities  are  comparable.  A  conservative 
upper  limit  on  the  velocity  of  the  scattering  material  is  1800  km  s-1,  based  on  the  lack  of  changes  in  the  shapes  of 
the  0.33  GHz  images.  In  the  0.33  GHz  field  of  view,  there  are  a  number  of  other  sources  that  might  also  be  heavily 
scattered,  which  would  suggest  that  there  are  large  changes  in  the  scattering  strength  on  lines  of  sight  separated  by 
a  degree  or  less.  Both  B 1849+005  and  PSR  B 1849+00  are  highly  scattered,  and  they  are  separated  by  only  13'.  If 
the  lines  of  sight  are  affected  by  the  same  “clump”  of  scattering  material,  it  must  be  at  least  2.3  kpc  distant. 
However,  a  detailed  attempt  to  account  for  the  scattering  observables  toward  these  sources  (angular  broadening  of 
the  extragalactic  source,  pulse  broadening  of  the  pulsar,  and  upper  limits  on  the  angular  broadening  of  the  pulsar) 
does  not  produce  a  self-consistent  set  of  parameters  for  a  clump  that  can  reproduce  all  three  measured  scattering 
observables.  A  clump  of  Ha  emission,  possibly  associated  with  the  H  n  region  G33.41 8— 0.004,  lies  between  these 
two  lines  of  sight,  but  it  seems  unable  to  account  for  all  of  the  required  excess  scattering. 

Subject  headings:  galaxies:  individual  (B 1849+005)  —  ISM:  structure  —  pulsars:  individual  (B 1849+00)  — 
radio  continuum:  general  —  scattering 


1.  INTRODUCTION 

Density  fluctuations  on  subparsec  scales  in  the  interstellar 
plasma  produce  refractive  index  fluctuations.  In  turn,  these 
refractive  index  fluctuations  scatter  radio  radiation  propagating 
through  the  plasma.  The  result  is  a  rich  set  of  observables 
including  angular  broadening  of  compact  sources,  pulsar  in¬ 
tensity  scintillations,  and  low-frequency  variability  of  extra¬ 
galactic  sources  (for  reviews  see  Rickett  1990;  Narayan  1992). 

The  density  fluctuations  responsible  for  interstellar  scatter¬ 
ing  appear  to  be  well  parameterized  by  a  power-law  power 
spectrum.  In  many  cases  the  spectral  index  of  the  power 
spectrum  is  near  11/3,  the  value  predicted  by  the  Kolmogorov 
theory  of  turbulence  for  neutral,  incompressible  unmagnetized 
fluids.  As  the  interstellar  plasma  is  compressible,  ionized,  and 
magnetized,  it  would  appear  to  violate  these  assumptions. 
There  has  been  recent  progress  in  understanding  how  a 
Kolmogorov-like  spectrum  could  arise  in  a  magnetized  plasma, 
however  (Lithwick  &  Goldreich  2003  and  references  within). 
Combined  with  the  large  Reynolds  number  of  the  interstellar 
plasma,  this  suggests  that  the  density  fluctuations  arise  from  a 
turbulent  process.  A  small  number  of  lines  of  sight  toward 
pulsars  are  notable  counterexamples,  having  spectral  indices 
larger  than  4  (Hewish  et  al.  1985;  Wolzszcan  &  Cordes  1987; 
Gupta  et  al.  1988;  Clegg  et  al.  1993)  or  an  index  approaching  4 
(Spangler  &  Gwinn  1990;  Moran  et  al.  1990;  Wilkinson  et  al. 
1994;  Molnar  et  al.  1995).  An  index  larger  than  4  arises  from  a 
random  distribution  of  density  fluctuations,  while  an  index 
approaching  4  is  interpreted  as  a  modification  of  the  power-law 
spectrum  by  an  inner  scale  indicative  of  a  dissipative  scale  for 


the  turbulence.  Thus,  determination  of  the  spectral  index  of 
the  density  fluctuation  spectrum  can  provide  clues  about  the 
mechanisms  by  which  the  density  fluctuations  originate. 

Wilkinson  et  al.  (1994)  showed  that  not  only  can  a  density 
spectral  index  be  obtained,  but  multifrequency  observations 
may  probe  the  organization  and  orientation  of  the  density 
fluctuations  responsible  for  scattering.  They  showed  that  the 
axial  ratio  and  position  angle  of  Cyg  X-3  are  frequency  de¬ 
pendent,  a  frequency  dependence  that  they  interpreted  in  terms 
of  an  increasingly  ordered  magnetic  field  on  smaller  and  smaller 
scales.  Trotter  et  al.  (1998)  showed  subsequently  that  the  axial 
ratio  and  orientation  of  the  scattering  disk  of  NGC  6334B  also 
exhibited  a  frequency  dependence.  Thus,  scattering  observa¬ 
tions  may  be  a  probe  of  the  magnetic  field  on  subparsec  scales. 

Conversely,  OH  masers  seen  toward  the  Galactic  center 
exhibit  some  of  the  most  extreme  anisotropic  images  seen 
anywhere  in  the  Galaxy  (van  Langevelde  et  al.  1992;  Frail  et  al. 
1994).  While  it  is  not  possible  to  measure  a  frequency  depen¬ 
dence  for  the  scattering  disks  of  the  OH  masers,  the  scattering 
disk  of  Sgr  A*  exhibits  no  change  in  its  axial  ratio  or  orien¬ 
tation  over  nearly  1.5  decades  in  frequency  (Lo  et  al.  1985; 
Jauncey  et  al.  1989;  Alberdi  et  al.  1993;  Yusef-Zadeh  et  al. 
1994;  Marcaide  et  al.  1999;  Reid  et  al.  1999).  It  is  not  yet  clear 
what  the  different  frequency  behavior  for  the  scattering  diam¬ 
eters  of  Cyg  X-3  and  NGC  6334B  versus  Sgr  A*  represents. 
Possibilities  might  include  a  difference  between  the  scattering 
media  in  the  Galactic  center  and  disk  or  a  much  stronger 
magnetic  field  in  the  Galactic  center. 

The  source  B 1849+005  (Galactic  coordinates  £  =  33?44, 
b  =  +0?21)  is  seen  along  the  third-most  highly  scattered  line  of 
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TABLE  1 

Observing  Log 


Frequency  Bandwidth  On-Source  Time  rms  Image  Noise 

(GHz)  Epoch  Array  (MHz)  Recorded  Polarization  (hr)  (mJy  beam-1) 


0.33 .  1987  Aug  20  VLA  (A-conf.)  3.0  R  0.28  80 

0.33 .  1998  Mar  19  VLA  (A-conf.)  3.0  R,  L  4.9  4.4 

1.5 .  1987  Aug  20  VLA  (A-conf.)  43.75  R  0.375  0.32 

2.3  .  2002  Aug  23  VLBA  (FD,  KP,  OV,  LA,  PT)  32  R,  L  1.2 

4.9 .  1987  Aug  20  VLA  (A-conf.)  50  R  0.017  1 

5.0 .  1998  Mar  29  VLBA+VLA1  56  L  1.4  0.39 

8.4  .  1998  Mar  29  VLBA+VLA1  64  R  0.89  0.51 

15 .  1998  Mar  29  VLBA+VLA1  64  L  1.4  0.84 

22 .  1998  Mar  29  VLBA+VLA1  64  L  1.4 


sight  known,  after  the  Galactic  center  (Lazio  &  Cordes  1998) 
and  NGC  6334B  (Trotter  et  al.  1998).  Spangler  et  al.  (1986) 
conclude  that  it  is  extragalactic  because  it  is  a  compact,  flat- 
spectrum  variable  source.  Further  evidence  that  B 1849+005  is 
extragalactic  is  found  by  comparing  its  H  i  absorption  spectrum 
to  that  of  the  pulsar  PSR  B 1849+00,  which  is  13'  from 
B 1849+005.  It  must  be  more  distant  than  PSR  B 1849+00  and 
is  probably  at  least  25  kpc  distant  (van  Gorkom  et  al.  1982; 
Dickey  et  al.  1983;  Clifton  et  al.  1988;  Frail  &  Weisberg  1990). 
PSR  B 1849+00  is  also  the  most  heavily  scattered  pulsar 
known,  with  a  pulse  broadening  time  of  0.2  seconds  at  1 .4  GHz. 
Thus,  a  comparison  of  the  lines  of  sight  to  B 1849+005  and 
PSR  B 1849+00  can  probe  changes  in  the  scattering  medium 
transverse  to  the  line  of  sight  on  scales  not  normally  accessible. 

This  paper  reports  multifrequency  observations  of  B 1849+ 
005,  spanning  nearly  two  decades  in  frequency,  that  are  used  to 
reassess  the  structure  and  nature  of  B 1849+005  and  constrain 
the  scattering  properties  along  this  line  of  sight.  In  §  2  I  de¬ 
scribe  the  observations  and  data  reduction,  in  §  3  I  use  these 
observations  to  extract  various  scattering  parameters  along  the 
line  of  sight  to  B 1849+005,  in  §  4  I  combine  my  observations 
with  those  from  the  literature  to  probe  the  scattering  environ¬ 
ment  around  B 1 849+005  and  particularly  toward  PSR  B 1 849+00, 
and  in  §  5  I  summarize  my  results. 

2.  OBSERVATIONS  AND  ANALYSIS 

The  scattering  disk  of  B 1849+005  is  large  enough  that  it  can 
be  resolved  by  a  range  of  frequencies  and  baselines  accessible 
with  the  Very  Large  Array  (VLA)  and  Very  Long  Baseline 
Array  (VLBA).  The  details  of  the  observations  and  data  re¬ 
duction  for  the  VLA  and  VLBA  observations  are  sufficiently 
different  that  I  discuss  them  separately.  Table  1  presents  the 
observing  log. 

2.1.  VLA  Observations  (0.33-4.9  GHz) 

2.1.1.  0.33  GHz 

Observations  were  obtained  at  0.33  GHz  with  the  A- 
configuration  of  the  VLA  on  1998  March  19.  The  sources  3C  48 
and  3C  286  were  used  in  calibrating  the  visibility  amplitudes, 
and  the  VLA  calibrator  J1 859+1259  was  used  in  calibrating 
the  visibility  phases.  Standard  calibration  and  radio  frequency 
interference  (RFI)  excision  procedures  were  used  within  the 
NRAO  Astronomical  Image  Processing  System  (AIPS).  Sev¬ 
eral  iterations  of  self-calibration  and  imaging  were  used  to 
approach  the  thermal  noise  level.  During  imaging,  a  poly¬ 
hedron  algorithm  (Cornwell  &  Perley  1992)  was  used  to  ac¬ 
count  for  the  noncoplanar  nature  of  the  VLA. 


At  0.33  GHz  the  field  of  view  is  large  (2?5  FWHM)  and 
contains  many  sources.  In  order  to  isolate  the  visibilities  from 
B 1849+005,  the  following  procedure  was  used.  Sources  within 
the  primary  beam  field  of  view  as  well  as  strong  sources  out¬ 
side  the  primary  beam  were  identified  in  the  NRAO  VLA  Sky 
Survey  (Condon  et  al.  1998).  Seventy-three  such  sources  were 
identified,  including  B 1849+005.  Following  Perley  (1989),  the 
Fourier  transform  of  all  of  the  sources  was  subtracted  from  the 
visibility  data,  a  modest  amount  of  additional  data  editing  was 
performed,  and  the  Fourier  transform  of  only  B 1849+005  was 
added  back  to  the  visibility  data. 

The  imaging  was  done  with  the  visibility  data  at  the  spectral 
resolution  at  which  it  was  acquired  (97.7  kHz  channels),  both  to 
avoid  bandwidth  smearing  and  to  make  RFI  easier  to  identify. 
After  removing  all  other  sources,  bandwidth  smearing  is  no  longer 
an  issue,  and  the  data  were  averaged  over  frequency.  Also,  in  order 
to  reduce  the  data  volume  before  further  processing,  the  visibili¬ 
ties  were  averaged  from  10  to  240  s.  This  time  and  frequency 
averaging  also  further  “dilutes”  the  impact  of  any  unsubtracted 
sources  or  errors  in  the  models  of  the  sources  that  were  subtracted. 

I  also  obtained  from  the  VLA  Archive  (program  AS305) 
A-configuration  VLA  observations  of  B 1849+005  acquired 
originally  on  1987  August  20  by  Fey  et  al.  (1991).  The  cali¬ 
bration  of  these  observations  was  similar  (but  with  the  VLA 
calibrator  B1 93 8  — 155  used  to  calibrate  the  phases),  but  given 
the  relatively  limited  amount  of  visibility  data  (15  minute  ob¬ 
servation),  only  B 1849+005  was  imaged,  rather  than  using 
targeted  facetting.  The  limited  visibility  data  also  implies  a 
higher  noise  level  and  limited  image  fidelity.  Figure  1  shows 
the  resulting  0.33  GHz  images  of  B1849+005. 

2.1.2.  1.5  GHz 

The  nominal  resolution  of  the  A-configuration  of  the  VLA 
at  1.5  GHz  is  approximately  L'5.  The  expected  scattering  di¬ 
ameter  of  B 1849+005  at  this  frequency  is  approximately  0'r3. 
“Superresolution”  imaging  can  provide  better  than  the  nom¬ 
inal  resolution  by  appropriate  weighting  of  the  visibility  data 
at  the  cost  of  a  modest  reduction  in  sensitivity  (Briggs  et  al. 
1999).  Given  the  strength  of  this  source  (~0.5  Jy  at  1.5  GHz), 
this  is  an  acceptable  trade-off. 

I  obtained  from  the  VLA  Archive  (program  AS305)  A- 
configuration  VLA  observations  of  B 1849+005  acquired 
originally  on  1987  August  20  by  Fey  et  al.  (1991)  with  the  aim 
of  imaging  the  source  at  superresolution.  The  source  3C  286 
was  used  to  calibrate  the  visibility  amplitude,  and  the  VLA 
calibrator  B 1821+107  was  used  in  calibrating  the  visibility 
phases.  Standard  calibration  procedures  were  used,  including 
a  number  of  iterations  of  phase  self-calibration. 
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Fig.  1. — B1849+005  at  0.33  GHz.  Left:  The  image  resulting  from  the  1987  August  20  observations.  The  rms  noise  level  is  80  mJy  beam-1,  and  the  contour  levels 
are  80  mJy  beam-1  x  —3,  3,  5,  7.07,  10,  ... .  The  beam  is  1577  x  571  and  is  shown  in  the  lower  left.  Right :  The  image  resulting  from  the  1998  March  19  observations. 
The  rms  noise  level  is  4.4  mJy  beam-1,  and  contour  levels  are  4.4  mJy  beam-1  x  —3,  3,  5,  7.07,  10, _ The  beam  is  774  x  479  and  is  shown  in  the  lower  left. 


As  was  done  in  imaging  the  0.33  GHz  observations,  a 
polyhedral  algorithm  was  used  to  account  for  the  non¬ 
coplanarity  of  the  VLA  and  image  the  entire  field  of  view  (30' 
FWHM).  In  practice,  the  number  of  other  sources  in  the  field 
was  sufficiently  small  and  their  flux  densities  sufficiently  low 
that  ignoring  the  VLA’s  noncoplanarity  would  have  made  little 
difference  in  the  final  image. 

Experiments  with  the  amount  of  superresolution  used 
showed  that  reasonable  choices  for  the  amount  of  weighting 
led  to  resolution  improvements  by  as  much  as  20%  (i.e.,  de¬ 
crease  in  beam  diameter).  Figure  2  shows  the  resulting  image 
of  B 1849+005  with  a  modest  amount  of  superresolution  (10% 
reduction  in  the  beam  diameter).  In  practice,  the  amount  of 
superresolution  used  made  little  difference  in  the  measured 
image  diameter,  but  the  measured  image  diameter  is  smaller 
than  the  beam  diameter.  Decreasing  the  beam  diameter  in¬ 
creases  the  ratio  of  the  image  diameter  to  the  beam  diameter 
and  therefore  increases  the  robustness  of  the  measured  image 
diameter,  particularly  if  the  image  diameter  is  smaller  than 
the  beam  diameter,  as  is  the  case  here. 

2.1.3.  4.9  GHz 

Fey  et  al.  (1991)  also  obtained  a  single  hour  angle  (1  minute 
duration)  on  B 1849+005  at  4.9  GHz  during  the  same  set  of 
observations  used  to  acquire  the  1.5  GHz  observations.  Based 
on  the  results  of  imaging  the  1.5  GHz  observations,  these 
4.9  GHz  observations  were  reanalyzed  also.  The  source 
3C  286  was  used  to  calibrate  the  visibility  amplitude,  and  the 
VLA  calibrator  B 1821+107  was  used  in  calibrating  the  visi¬ 
bility  phases.  Standard  calibration  procedures  were  used,  in¬ 
cluding  a  small  number  of  iterations  of  phase  self-calibration. 
Figure  3  shows  the  resulting  image. 

2.2.  VLB  A  Observations  (2.3-22  GHz) 

Observations  at  5,  8.4,  15,  and  22  GHz  were  obtained  with 
the  VLB  A  and  one  antenna  of  the  VLA  on  1998  March  29. 
Observations  at  2.3  GHz  were  obtained  with  only  the  inner 
antennas  of  the  VLB  A  on  2002  August  23.  No  fringes  were 


detected  at  22  GHz.  The  lack  of  a  detection  at  22  GHz  is 
consistent  with  the  lower  sensitivity  of  the  VLBA  at  this 
frequency  relative  to  the  other  frequencies  observed.  At  the 
other  frequencies,  the  correlated  visibilities  were  fringe  fit  and 
amplitude  calibrated  in  the  standard  fashion  within  AIPS. 
Fringe-fitting  intervals  were  2.5-3  minutes.  The  calibrated  vis¬ 
ibilities  were  then  exported  to  the  Caltech  Difmap  program 
for  hybrid  mapping. 


Fig.  2. — B 1849+005  at  1.5  GHz,  imaged  using  a  modest  amount  of  super¬ 
resolution  to  obtain  a  10%  reduction  in  the  nominal  beam  diameter.  The  rms 
noise  level  is  0.32  mJy  beam-1,  and  contour  levels  are  0.32  mJy  beam-1  x  — 3, 
5,  7.07,  10,  14.1,  20, _ The  beam  is  173  x  173  and  is  shown  in  the  lower  left. 
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Fig.  3. — B 1849+005  at  4.9  GHz,  as  observed  by  the  VLA.  The  rms  noise 
level  is  1  mJy  beam-1,  and  contour  levels  are  1  mJy  beam-1  x  —3,  5,  7.07, 
10,  14.1,  ....  The  beam  is  0"5  x  0"4  and  is  shown  in  the  lower  left. 


For  the  5,  8.4,  and  15  GHz  observations,  after  obtaining  an 
image  with  a  noise  near  but  still  well  above  the  thermal  limit, 
I  subtracted  the  Fourier  transform  of  B 1849+005  from  the 
visibility  data,  performed  a  modest  amount  of  additional  data 
editing,  and  then  added  the  Fourier  transform  of  B 1849+005 
back  to  the  visibility  data  (Perley  1989).  Figures  4-6  show  the 
resulting  images  of  B 1849+005  between  5  and  15  GHz. 

At  2.3  GHz,  a  significant  detection  of  B 1849+005  was 
obtained  only  on  the  shortest  baseline,  between  the  LA  and  PT 
antennas.  As  shown  below,  this  source  has  a  fairly  simple 
structure  on  milliarcsecond  scales.  In  order  to  estimate  a  di¬ 
ameter,  the  correlated  visibility  amplitude  on  the  LA-PT 
baseline  was  compared  to  total  flux  density  of  the  source.  Furst 
et  al.  (1990)  measure  the  flux  density  of  the  source  at  2.7  GHz 
to  be  0.77  Jy;  flux  densities  obtained  by  interpolating  its 
spectrum  between  1.5  and  5  GHz  (Fey  et  al.  1991;  this  work) 
range  from  0.6  to  0.8  Jy.  Because  the  measurement  by  Furst 
et  al.  (1990)  utilized  a  single-dish  telescope,  I  adopt  0.77  Jy 
as  its  flux  density  at  2.3  GHz.  Assuming  the  source  to  be  a 
circular  Gaussian,  its  diameter  is  then  110  mas  (Table  2). 

For  further  analysis,  the  data  were  averaged  in  time.  The 
averaging  times  were  30  s  at  2.3  GHz,  60  s  at  5  GHz,  30  s 
at  8.4  GHz,  and  15  s  at  15  GHz. 

3.  INTERSTELLAR  SCATTERING  TOWARD  B 1849+005 
3.1.  Relevant  Scattering  Formulae 

The  interferometric  visibility  of  an  infinitely  distant  point 
source  seen  through  a  medium  filled  with  random  density 
fluctuations  is 

F(Z>)  =  e_A+)/2  (1) 

for  an  interferometer  with  baseline  length  b.  For  the  density 
power  spectrum  relevant  in  the  interstellar  medium  (Armstrong 


Fig.  4. — Central  component  of  B 1849+005  at  5  GHz,  as  observed  by  the 
VLBA.  The  rms  noise  level  is  0.39  mJy  beam-1,  and  contour  levels  are 
0.39  mJy  beam-1  x  —2,  3,  5,  7.07,  10,  ....  The  beam  is  14  x  7.6  mas  and  is 
shown  in  the  lower  left. 

et  al.  1995),  the  form  of  the  phase  structure  function  depends 
upon  the  value  of  the  importance  of  l\  relative  to  the  typical 
baseline  length  b.  Anticipating  later  results,  I  consider  two 
possibilities  for  D^b)  (e.g.,  Coles  et  al.  1987;  Cordes  &  Lazio 
1991;  Lambert  &  Rickett  1999) 

D0(b)  =  8tt r2eX2  SM 

f(a)ba~2  b»lu 

f(a)T(a/2)  l^F^a /2~  1;  -(b/hf)  -  1]  b<lx. 

(2) 


Fig.  5. — Central  component  of  B 1849+005  at  8.4  GHz.  The  rms  noise  level 
is  0.51  mJy  beam-1,  and  contour  levels  are  0.51  mJy  beam-1  x  —3,  3,  5,  7.07, 
10, _ The  beam  is  4.2  x  2.8  mas  and  is  shown  in  the  lower  left. 
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Fig.  6. — Central  component  of  B 1849+005  at  15  GHz.  The  rms  noise  level 
is  0.84  mJy  beam-1,  and  contour  levels  are  0.84  mJy  beam-1  x  —2,  3,  5,  7.07, 
10, _ The  beam  is  2.0  x  0.88  mas  and  is  shown  in  the  lower  left. 


Here  re  is  the  classical  electron  radius,  X  is  the  free-space 
observing  wavelength, /(ct)  is  a  slowly  varying  function  of  a 
with  a  value  near  unity,  \F\(a\b\x)  is  the  confluent  hyper¬ 
geometric  function,  and  the  scattering  measure  SM  is  pro¬ 
portional  to  the  line-of-sight  integrated  rms  electron  density. 
(Minor  modifications  are  required  if  the  source  is  not  infin¬ 
itely  distant,  as  assumed  here.)  For  reference,  in  the  local  in¬ 
terstellar  medium,  SM  «  10-3  5  kpc  m-20/3  for  a  1  kpc  path 
length. 

The  resulting  scattering  diameter,  for  a  distant  source  and 
assuming  a  =  11/3  (Rickett  1990;  Cordes  &  Lazio  1991),  is 


0d 


128  mas  SM0  6 
71  mas  SM1/2 


b  »/i, 

^GHzf 

h  \ 

-1/6 

b<£lXl 

100  km/ 

(3) 


where  z^ghz  is  the  observing  frequency  in  GHz  and  SM  is  in 
its  canonical  units  of  kpc  m_20/3. 

3.2.  Angular  Broadening  and  the  Apparent  Structure 
of  B 1849+ 005 

On  the  basis  of  its  spectrum,  Fey  et  al.  (1991)  suggested  that 
the  structure  of  B 1849+005  may  be  more  complex  than  that  of 
a  single  compact  component.  Figures  2  and  3  demonstrate  that 
this  is  the  case.  Although  a  strong  central  component  does 
dominate  the  structure  of  the  source,  there  are  weaker  com¬ 
ponents  both  to  the  southwest  and  northeast.  The  components 
to  the  northeast  are  at  best  only  faintly  visible  in  the  VLA 
observations  at  4.9  GHz  (Fig.  3),  though  their  absence  prob¬ 
ably  results  from  the  limited  hour-angle  coverage  and  surface 
brightness  sensitivity  of  the  A-configuration  of  the  VLA  at  this 
frequency.  The  central  component  may  show  a  hint  of  a  jet 
emerging  to  the  southeast  and  bending  toward  the  southwest. 
The  existence  of  a  jet  must  be  considered  tentative,  however, 
given  the  extremely  limited  hour-angle  coverage  (only  1  minute 
observation).  As  the  0.33  GHz  image  is  elongated  in  approx¬ 
imately  the  same  orientation  as  that  of  the  central  and  south¬ 
west  components,  I  have  also  reimaged  the  (1998  March  19) 
0.33  GHz  observations  utilizing  superresolution.  There  is  no 
indication  of  this  second  component,  though  the  orientation 
of  the  beam  is  similar  to  that  of  the  central  and  southwest 
components. 

The  expressions  of  §  3.1  assume  an  infinitely  distant  point 
source.  To  what  extent  can  they  be  applied  to  B 1849+005, 
given  its  complex  structure?  Figure  7  shows  the  diameter  of  the 
central  component  of  B 1 849+005  as  a  function  of  observing 
frequency.  Although  Figures  1-6  show  images  of  B 1849+005, 
diameters  measured  from  these  images  were  not  used  in  con¬ 
structing  Figure  7.  Rather,  Gaussian  models  were  fitted  to  the 
visibility  data,  except  in  the  case  of  the  2.3  GHz  observations 
for  which  a  circular  Gaussian  was  assumed  (§  2).  Also  plotted 
in  Figure  7  are  measurements  at  various  other  frequencies.  At 
these  other  frequencies  (0.25  and  0.408  GHz)  the  visibility  data 
were  not  available  for  analysis  (as  below),  but  angular  diam¬ 
eters  could  be  determined.  For  completeness,  Table  2  sum¬ 
marizes  the  various  fits  for  the  diameter  of  the  source. 

Figure  7  also  shows  three  curves,  0  oc  v~2-2 ,  0  oc  v~2 ,  and 
0  oc  v~XA .  The  first  two  are  relevant  for  strongly  scattered 
sources,  with  the  first  being  appropriate  if  the  electron  density 


TABLE  2 

Gaussian  Model  Fits  to  the  Central  Component  of  B 1849+005 


Frequency  S  a  b  § 

(GHz)  (Jy)  (mas)  (mas)  (deg)  References 


0.25 .  1.624  +  0.009  11,400  +  70  8430  +  60  44.2  +  0.7  1 

0.33 .  2.1  ±0.2  6300  ±  1000  2080  ±100  18  ±7  2 

0.33 .  2.12  ±0.02  7050  +  350  2560+  100  25.8  ±0.7  2 

0.41 .  1.7  ±0.2  >600  ...  ...  3 

1.5a .  0.63  ±0.06  440  +  6  190  +  6  161  ±  1  2 

2.3  .  0.77  +  0.1  110  +  5  110  +  5  ...  2 

5.0 .  0.739  ±0.009  19.34+  1  17.32+  1  27.0  ±0.6  2 

8.4  .  0.634  ±0.006  6.71  ±0.3  4.91  ±0.2  -33.2  ±0.3  2 

15 .  0.831  ±0.008  2.017  +  0.1  1.623  +  0.08  -21.1  ±0.5  2 


Notes. — The  visibility  data  are  fit  with  a  single  Gaussian  model  having  a  flux  density  S,  major  axis  a,  minor  axis  b,  and 
position  angle  4>.  In  some  cases  a  circular  Gaussian  was  fit;  no  position  angle  is  reported  for  these  cases.  The  first  entry 
for  0.33  GHz  is  for  the  1987  August  20  observations;  the  second  is  for  the  1998  March  19  observations.  The  stated  uncertainties 
are  the  formal  statistical  uncertainties  for  the  fits;  the  actual  uncertainties  due  to  systematic  effects  are  almost  certainly  larger. 
a  Only  the  central  component  was  fitted. 

References. — (1)  A.  R  Rao  (2002,  private  communication);  (2)  this  work;  (3)  Dennison  et  al.  (1984). 
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Fig.  7. — Angular  diameter  of  B 1849+005  as  a  function  of  observing  fre¬ 
quency.  The  angular  diameter  was  obtained  by  fitting  a  Gaussian  either  to  the 
image  or  to  the  visibility;  if  an  elliptical  Gaussian  was  used  in  the  fit,  shown  is 
the  geometric  mean  of  the  major  and  minor  axes.  At  most  frequencies,  the 
formal  uncertainties  are  comparable  to  the  size  of  the  points.  The  datum 
at  0.25  GHz  is  from  P.  Rao  (2001,  private  communication)  and  that  at  0.408 
from  Dennison  et  al.  (1984).  For  reference,  three  curves  are  also  plotted,  one 
showing  a  v~2  dependence  for  the  diameter  ( dotted  line),  one  showing  a  is~2'2 
dependence  ( dashed  line),  and  one  showing  a  v~lA  dependence  ( dot-dashed 
line)  that  would  be  expected  if  the  angular  broadening  for  B 1849+005  is 
anomalous  in  a  manner  similar  to  that  for  PSR  B 1849+00  (§  4.2).  All  three 
curves  are  constrained  to  go  through  the  diameter  measured  at  5  GHz. 

spectrum  is  a  power  law  with  a  Kolmogorov  spectral  index, 
a  =  11/3,  while  the  second  is  appropriate  if  the  range  of 
observing  baselines  is  comparable  to  l\.  The  third  curve  shows 
the  frequency  dependence  for  the  angular  broadening  if  scat¬ 
tering  toward  B 1849+005  is  anomalous  in  a  manner  similar  to 
that  seen  for  the  pulsar  PSR  B 1849+00  (§  4.2).  The  relatively 
good  agreement  between  the  angular  diameters  and  the  first 
two  curves  is  an  indication  that  scattering  dominates  the  ap¬ 
parent  structure  of  the  central  component  at  least  up  to  15  GHz. 
Interpolating  the  angular  diameter  to  the  fiducial  frequency 
of  1  GHz,  the  scattering  diameter  of  B 1849+005  is  0"6.  Based 
on  their  more-limited  set  of  observations,  Fey  et  al.  (1991) 
reached  a  similar  conclusion:  B 1849+005  is  heavily  scattered 
but  exhibits  a  more  complex  structure  than  simply  a  single 
compact  component. 

In  addition,  the  angular  diameter  of  the  southwest  component 
was  measured,  using  the  images  shown  in  Figures  2  and  3.  Its 
diameter  is  approximately  constant  between  1.5  and  4.9  GHz, 
indicating  that  its  diameter  is  much  larger  than  the  scattering 
diameter.  This  is  consistent  with  its  somewhat  diffuse  appear¬ 
ance  at  4.9  GHz. 

3.3.  The  Electron  Density  Spectrum 

A  more  profitable  analysis  of  the  observations  is  to  use  the 
structure  function  (eqs.  [1  ]— [2])  to  extract  parameters  of  the 
electron  density  spectrum.  My  analysis  follows  closely  that  of 
Spangler  &  Cordes  (1998).  Figure  8  shows  an  example  of 
a  structure  function  constructed  from  the  8.4  GHz  visibility 
data. 

Clearly  apparent  in  Figure  8  is  a  nonrandom  scatter  in  the 
structure  function,  e.g.,  around  10  M2.  Moreover,  Figures  1 
and  4-6  show  that  the  scattered  image  has  (varying)  degrees 
of  anisotropy.  This  anisotropy  arises  presumably  because  the 
density  fluctuations  responsible  for  the  scattering  are  them¬ 
selves  anisotropic. 


T - 1 - 1  I  | 
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Fig.  8. — Structure  function  for  the  8.4  GHz  visibility  data  using  the  iso¬ 
tropic  baseline.  A  small  number  of  visibilities  below  2  M2  are  not  shown. 

The  anisotropy  of  the  image  and  visibility  data  (and  of 
the  underlying  density  fluctuations)  are  described  by  an 
orientation  and  aspect  ratio  r).  In  fitting  the  visibility  data  to 
equation  (1),  the  baseline  b  is  replaced  by  the  “rotundate” 
baseline 

b2r  =  ( u  cos  +  v  sin  \f)2  +  rj2(—  u  sin  +  v  cos  4>)2 ,  (4) 

where  u  and  v  are  the  traditional  interferometer  baseline 
coordinates  measured  in  units  of  the  observing  wavelength. 
The  rotundate  baseline  is  simply  the  baseline  measured  in  an 
anisotropic  coordinate  system  that  is  rotated  relative  to  that  of 
the  interferometer’s  coordinate  system.  In  this  rotundate  co¬ 
ordinate  system,  4/  is  measured  counterclockwise  from  the 
u  axis,  in  contrast  to  the  position  image  in  an  image,  which  is 
measured  counterclockwise  from  north.  Thus,  there  is  a  90° 
offset  from  an  image  position  angle  and  \F. 

The  nonrandom  scatter  in  Figure  8  arises  because  of  the  use 
of  the  isotropic  baseline  (77  —  1  and  T  =  0°).  Figure  9  illus¬ 
trates  the  use  of  the  rotundate  baseline  in  plotting  the  structure 
function  (using  the  best-fitting  power-law  density  spectrum 
model  found  below).  Clearly  apparent  is  the  dramatic  reduc¬ 
tion  in  the  scatter  of  the  data;  similar  improvements  are  seen  at 
the  other  frequencies.  Figure  9  also  demonstrates  that  the 
structure  of  B 1 849+005  does  not  have  an  adverse  impact  on 
efforts  to  extract  information  about  the  scattering.  A  single 
source  component  does  an  excellent  job  of  fitting  the  visibil¬ 
ities.  As  might  be  expected  from  Figures  2  and  3,  the  VLBI 
observations  are  insensitive  to  the  relatively  low  brightness 
temperature  components  to  the  southwest  and  northeast. 

One  important  aspect  of  the  use  of  the  rotundate  baseline  is 
that  the  structure  functions  in  equation  (2)  have  been  derived 
assuming  an  isotropic  density  spectrum.  This  ex  post  facto 
replacement  of  the  rotundate  baseline  for  the  baseline  violates 
that  assumption.  This  replacement  should  be  approximately 
valid  provided  that  77  is  not  much  larger  than  unity.  Moreover, 
Chandran  &  Backer  (2002)  have  considered  the  propagation 
of  radio  waves  through  an  anisotropic  medium.  They  find  that 
expressions  like  those  of  equation  (2)  continue  to  hold;  in 
particular,  the  power-law  dependences  remain  unchanged. 

Using  rotundate  baselines  and  the  structure  functions  of 
equation  (2),  there  are  four  or  five  parameters  for  which  to  fit, 
depending  upon  the  relative  importance  of  the  inner  scale:  a, 
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Fig.  9. — Structure  function  for  the  8.4  GHz  visibility  data  using  the 
rotundate  baseline.  The  solid  line  shows  the  best-fitting  model  for  a  power-law 
density  spectrum  (viz.,  Table  3).  A  small  number  of  visibilities  below  2  M2 
are  not  shown. 

SM,  r],  and  l\.  I  used  a  genetic  algorithm  (Lazio  1997) 
to  find  the  best-fitting  parameters  in  a  (unweighted)  minimum 
X2  sense  at  each  frequency. 

I  have  fit  the  structure  functions  (as  in  Fig.  8)  for  the  VLB  I 
observations  (2.3-15  GHz)  and  the  0.33  GHz  VLA  observa¬ 
tions  of  1998  March  19.  The  visibility  data  from  the  1.5  or 
4.9  GHz  observations  from  the  VLA  were  not  used,  because 
it  is  clear  that  the  source  structure  is  complex,  and  at  4.9  GHz, 


the  expected  scattering  diameter  is  much  smaller  than  what 
can  be  resolved  by  the  VLA.  In  addition,  the  0.33  GHz  visi¬ 
bility  data  from  the  1987  August  20  VLA  observations  were 
not  used  because  of  their  limited  quantity  and  restricted  u-v 
coverage. 

Not  all  of  the  data  were  used  in  the  fitting.  At  long  baselines, 
the  structure  functions  saturate,  while  at  short  baselines  they 
reach  noise  plateaus  (viz.,  eqs.  [1]  and  [2]).  The  saturation  at 
long  baselines  is  visible  near  30  M2  in  Figures  8  and  9;  the  noise 
plateau  at  short  baselines  occurs  below  2  M2  (for  the  8.4  GHz 
data).  Consequently,  limits  were  set  at  both  short  and  long 
baselines.  After  excluding  noisy  data,  the  total  number  of  data 
available  to  fit  ranged  from  several  hundred  to  several  thousand. 

Figures  10  and  11  show  the  best-fitting  parameters  and  their 
confidence  ranges  as  a  function  of  frequency  for  both  a  power- 
law  density  spectrum  (b  l\)  and  a  density  spectrum  with  an 
inner  scale  respectively;  Tables  3  and  4  tabulate  the 

best-fitting  model  parameters.  Figures  10  and  11  also  indicate 
the  (weighted)  mean  values,  over  all  frequencies,  for  the 
various  parameters.  A  x2  statistic  was  used  to  test  whether  the 
best-fit  values  are  consistent  with  a  single  mean  value,  i.e.,  no 
frequency  dependence. 

The  confidence  ranges  for  the  parameters  have  been  esti¬ 
mated  both  by  a  local  grid  search  and  a  nonlinear  least-squares 
minimization  around  the  minimum  %2.  In  the  former  method, 
the  confidence  range  for  a  parameter  was  determined  by 
varying  that  parameter,  while  holding  all  others  fixed  at  their 
best-fitting  values,  until  the  %2  approximately  doubled  from  its 
minimum  value.  In  the  latter  method  a  Levenberg-Marquardt 
minimization  routine  was  used  to  determine  the  covariance 


Fig.  10. — Best-fit  parameters  as  a  function  of  frequency  for  a  power-law  density  spectrum.  The  confidence  regions  shown  are  the  amount  by  which  the  parameter 
must  change  in  order  to  cause  the  x2  to  increase  by  unity  from  its  minimum  value.  The  horizontal  dotted  line  shows  the  weighted  mean  value  for  the  parameter. 
Because  B 1849+005  was  detected  on  only  a  single  baseline  at  2.3  GHz  {asterisk),  the  anisotropy  g  and  its  orientation  were  constrained  to  be  1  and  0  (cos  =  1), 
respectively.  At  1.5  GHz,  a  structure  function  analysis  was  not  conducted;  the  values  shown  are  derived  from  analysis  of  the  image  {triangle). 
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Fig.  11. — Same  as  Fig.  10,  but  for  best-fit  parameters  as  a  function  of  frequency  for  a  density  spectrum  with  an  inner  scale. 


matrix  of  the  parameters.  The  Levenberg-Marquardt  minimi¬ 
zation  yields  confidence  ranges  much  smaller  than  the  grid 
search.  We  have  adopted  the  confidence  ranges  for  the  grid 
search;  adopting  smaller  confidence  ranges  would  make  our 
following  conclusions  more  robust. 

One  aspect  of  Figures  10  and  11  warrants  comment.  The 
various  parameters  (a,  SM,  77,  and  l\)  are  plotted  as 
functions  of  frequency.  These  parameters  were  introduced  as 
frequency-independent  quantities  (eq.  [2]).  For  instance,  SM 
is  proportional  to  the  rms  electron  density.  The  apparent  fre¬ 
quency  dependence  arises  from  the  use  of  interferometers 
whose  baseline  lengths  were  not  scaled  with  frequency.  For  an 
interferometer  with  a  typical  (and  fixed)  baseline  length  b ,  it 
will  resolve  the  scattering  disk  only  if  the  scattering  disk  is 
comparable  in  size  to  the  typical  angle  measured  by  the  in¬ 
terferometer  kjb.  In  turn,  if  the  scattering  medium  is  at  a 

TABLE  3 

Power-Law  Density  Spectrum  Model  Fitting 


Frequency  ^  log  SM 

(GHz)  a  7 \  (deg)  (kpc  m_20/3) 


0.33 .  3.68^0  qi  2.29^-f  158+44  0.47+g^ 

1.5a .  ...  '  2.3’  161  0.98 

2.3b .  3.60  ...  ...  1.3 

5.0 .  3.94+QQ2  1.1 8+2:  is  134^  0.60^0^5 

8.4 .  3.83^:?2  1.49+™*  119+J3  0.79+™* 

15 .  3.91+q’q3  1.29  114+2b  0.73+™2 

Mean .  3.92  1.30  141  0.68 


a  At  this  frequency,  the  scattering  parameters  were  derived  from  the  image, 
not  from  a  structure  function  analysis. 

b  At  this  frequency,  the  scattering  diameter  was  assumed  to  be  circular. 


distance  D,  the  scattering  disk  will  be  produced  by  density 
fluctuations  on  scales  of  l'  ~  D(A/b).  This  length  scale  is 
manifestly  frequency  dependent.  Thus,  an  apparent  frequency 
dependence  in  a,  SM,  77,  T',  or  l\  would  reflect  different  pa¬ 
rameter  values  being  obtained  on  different  length  scales  within 
the  spectrum  of  scattering  irregularities.  For  brevity,  in  the 
context  of  Figures  10  and  11,  I  shall  continue  to  speak  of  a 
frequency  dependence  for  a,  SM,  77,  4/,  and  l\  with  the  un¬ 
derstanding  that  any  frequency  dependence  arises  from  the 
different  length  scales  probed  in  the  scattering  medium. 

3.3.1.  Image  Anisotropy 

Regardless  of  whether  a  power-law  or  inner-scale  density 
spectrum  is  fitted  to  the  structure  functions,  the  result  confirms 
what  the  casual  inspection  of  Figures  1-6  suggest:  the  scat¬ 
tering  toward  B 1849+005  is  decidedly  anisotropic. 

At  0.33  GHz,  the  anisotropy  appears  large,  77  ^  2.3,  but  the 
weaker  components,  especially  the  one  to  the  southeast,  may 


TABLE  4 

Power-Law  with  Inner-Scale  Density  Spectrum  Model  Fitting 


Frequency 

$ 

log  SM 

h 

(GHz) 

a 

1 

(deg) 

(kpc  m_20/3) 

(km) 

0.33 . 

3  70+™1 
/u-0.10 

2.00+?;§5 

160^7 

0-921“;“ 

1901^° 

2.3a . 

3.60 

1.301“;! 

62-u 

5.0 . 

8.4 . 

15 . 

3-80+8;?*, 

/U— 0.06 
3.80±°;°J 

1  20+010 
I.ZU— 020 

1  30+030 
1-:>u-0.20 

1  20+™4 
I.ZU— 020 

142+IS 

1251“ 

1  ?  1  +32 

1  14+0-21 
1-i4-0.42 

1  O3+0-25 

1-UJ-0.25 

1  16+004 
1>iO-0.58 

500!?00 

4701“ 

490!!q 

Mean . 

3.78 

1.23 

148 

1.14 

470 

a  At  this  frequency,  the  scattering  diameter  was  assumed  to  be  circular. 
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be  contaminating  the  analysis.  The  beam  diameter  at  0.33  GHz 
is  comparable  to  the  separation  between  the  components,  and 
the  orientation  of  the  source  at  1.5  GHz  is  similar  to  that  of 
both  the  source  and  the  position  angle  of  the  beam  at  0.33  GHz. 
The  visibilities  at  0.33  GHz  appear  consistent  with  only  a 
single  component  being  present  (§  3.2  and  Fig.  7),  but  the 
presence  of  the  other  components  cannot  be  excluded.  Indeed, 
convolving  the  1.5  GHz  image,  Figure  2,  to  a  resolution  com¬ 
parable  to  that  of  the  0.33  GHz  images,  Figure  1,  produces  an 
image  similar  to  the  0.33  GHz  images. 

Even  if  the  analysis  at  0.33  GHz  is  disregarded  (and  at  1.5 
and  2.3  GHz,  for  which  less  comprehensive  analyses  could  be 
performed),  the  higher- frequency  data  show  a  clear  anisotropy. 
The  fits  at  5.0  and  1 5  GHz  are  just  marginally  consistent  with  an 
isotropic  image,  while  the  fit  at  8.4  GHz  excludes  an  isotropic 
image.  The  mean  value  for  the  image  elongation  is  77  =  1.2. 
This  value  is  comparable  to  that  measured  for  Cyg  X-3  (~1.3, 
Wilkinson  et  al.  1994;  Molnar  et  al.  1995)  and  NGC  6334B 
(1.2- 1.5;  Trotter  et  al.  1998),  within  the  range  of  anisotropies 
measured  for  the  extragalactic  sources  seen  through  the  Cygnus 
region  (1.1-1. 8;  Spangler  &  Cordes  1988,  1998;  Desai  &  Fey 
2001),  and  somewhat  smaller  than  the  value  measured  for  the 
OH  masers  toward  W49N  (2-3;  Desai  et  al.  1994)  and  the 
Galactic  center  («2.5;  Frail  et  al.  1994).  Romani  et  al.  (1986) 
have  predicted  the  anisotropy  that  would  be  induced,  even  if  the 
small-scale  (diffractive)  density  fluctuations  responsible  for 
angular  broadening  are  isotropic,  by  refractive  effects  from 
large-scale  density  fluctuations.  For  a  thin  scattering  screen 
with  a  Kolmogorov  density  spectrum,  the  expected  anisotropy 
is  (77)  <  1.05.  The  scattering  toward  B 1849+005  is  both  strong 
and  anisotropic,  and  the  anisotropy  must  result  from  the  density 
fluctuations  responsible  for  the  scattering. 

Whether  the  position  angle  4/  of  the  anisotropy  changes 
with  frequency,  as  Wilkinson  et  al.  (1994)  and  Trotter  et  al. 
(1998)  found  for  Cyg  X-3  and  NGC  6334B,  respectively,  is 
less  certain.  In  both  Figures  10  and  11,  an  apparent  trend  of  a 
decreasing  position  angle  (or  increasing  cos  \h)  is  seen.  This  is 
the  case  regardless  of  whether  the  0.33  GHz  datum  is  included 
or  not.  However,  as  both  figures  show,  the  uncertainties  are 
also  consistent  with  no  change  in  position  angle. 

Using  only  the  higher- frequency  data,  the  scattering  diameter 
of  B 1849+005  ranges  from  19  mas  at  5.0  GHz  to  2  mas  at 
15  GHz.  If  the  material  responsible  for  the  scattering  is  at  a 
distance  of  5  kpc  (§  4.2),  the  size  of  the  region  on  which  the 
anisotropic  density  fluctuations  occur  must  range  from  1 .5  x  1014- 
1.4  xlO15  cm  (10-95  AU).  These  scales  are  comparable  to 
those  seen  for  NGC  6334B,  though  this  is  due  partially  to  a 
selection  effect  in  that  the  same  telescope  was  used  in  both  this 
work  and  that  of  Trotter  et  al.  (1998),  so  the  range  of  accessible 
spatial  scales  is  the  same. 

Within  the  context  of  the  models  for  magnetohydrodynamic 
turbulence  developed  by  Goldreich  and  collaborators  (e.g., 
Fithwick  &  Goldreich  2001),  these  length  scales  would  rep¬ 
resent  the  lower  limit  to  their  “MHD  scale.”  On  this  length 
scale,  the  kinetic  energy  density  (oc  v2)  of  the  turbulent  eddies 
has  diminished  to  the  point  that  the  magnetic  energy  density 
(oc  B 2)  is  comparable  to  it.  Thus,  this  length  scale  marks  where 
the  turbulence  makes  a  transition  from  isotropic  on  large  scale 
(comparable  to  the  outer  scale)  to  increasingly  anisotropic  on 
small  scales. 

3.3.2.  Inner  Scale 

Comparison  of  Figures  10  and  11  shows  that  the  density 
spectral  index  a  appears  to  change  between  0.33  and  5.0  GHz 


for  the  power-law  spectrum.  At  0.33  GHz,  the  fit  for  a  power- 
law  spectrum  finds  a  =  3.68,  while  above  5.0  GHz,  a  >  3.8. 
In  contrast,  the  fit  for  an  inner-scale  spectrum  reveals  a  much 
more  consistent  set  of  values  for  a,  with  3.7  <  a  <  3.8. 
Moreover,  the  fits  for  the  inner  scale  itself  find  consistent  values 
of  a  few  hundred  kilometers  (though  with  large  uncertainties). 

I  regard  this  as  a  detection  of  the  effects  of  the  inner  scale  on 
scales  of  a  few  hundred  kilometers.  The  maximum  baselines  of 
the  VFA,  used  for  the  0.33  GHz  observation,  are  35  km.  In 
contrast,  the  VFBA  baselines,  used  for  the  higher-frequency 
observations,  span  the  range  from  a  few  hundred  to  a  few 
thousand  kilometers.  Thus  the  “break”  between  0.33  and 
5.0  GHz  in  the  fits  for  a  in  the  power-law  density  spectrum 
(Fig.  10)  reflects  the  range  of  baselines  used  in  the  observations. 

This  value  for  the  inner  scale,  a  few  hundred  kilometers,  is 
comparable  to  that  found  toward  other  sources  (Spangler  & 
Gwinn  1990;  Moran  et  al.  1990;  Wilkinson  et  al.  1994;  Molnar 
et  al.  1995). 

3.4.  Secular  Change  of  the  Scattered  Image 

If  the  plasma  responsible  for  the  scattering  has  a  bulk 
(transverse)  velocity  v  (and  it  remains  relatively  unchanged  as 
it  moves  across  the  line  of  sight,  “frozen-flow”  approxima¬ 
tion),  changes  in  the  image  appearance  should  appear  on 
timescales  comparable  to  the  refractive  timescale  tr  =  lr/v. 
The  two  observations  at  0.33  GHz  (Fig.  1)  provide  a  10.6  yr 
baseline  over  which  to  search  for  secular  changes  in  the  image. 

Table  2  shows  that  both  the  fitted  diameters  and  position  angles 
for  the  source  do  not  change  appreciably.  From  both  sets  of 
observations  —  4"  at  0.33  GHz.  The  refractive  scale  is  there¬ 
fore  6  x  1016  cm  Dk pC,  where  Z>kpc  is  the  distance  to  the  screen  in 
kiloparsecs.  There  is  little  change  in  the  image  over  the  10.6  yr 
interval.  The  velocity  of  the  scattering  material  is  therefore 
v  <  1800  km  s-1  DkpC.  (Strictly,  this  velocity  is  a  combination 
of  the  velocities  of  the  source,  medium,  and  observer.  Here  the 
velocity  of  the  source  is  assumed  to  be  effectively  zero.) 

This  upper  limit  is  conservative.  If  the  actual  velocity  were 
within  a  factor  of  two,  say,  of  this  upper  limit,  some  modest 
changes  in  at  least  the  orientation  should  start  to  become  apparent. 

The  source  B 1849+005  is  only  14'  from  the  supernova 
remnant  (SNR)  G33. 6+0.1.  Spangler  et  al.  (1986)  observed 
B 1849+005  with  the  aim  of  assessing  to  what  extent  the  SNR 
could  contribute  to  the  scattering  of  B 1849+005.  While  they 
inferred  correctly  that  B 1849+005  is  heavily  scattered,  they 
were  reluctant  to  attribute  all  of  the  scattering  to  G3 3. 6+0.1 
given  the  low-latitude  line  of  sight  to  B 1849+005.  Unfortu¬ 
nately,  this  upper  limit  on  the  velocity  of  the  scattering  material 
sheds  no  light  on  whether  material  from  the  SNR  is  responsible 
for  the  scattering.  The  large  scattering  diameter  at  0.33  GHz 
means  that  relatively  large  length  scales  are  probed  within  the 
scattering  material  and  that  correspondingly  long  time  intervals 
are  needed  to  place  stringent  constraints  on  the  velocity. 

The  higher-frequency  VFBI  observations  may  prove  useful 
in  constraining  the  velocity  of  the  scattering  material  in  the 
future.  For  instance,  repeating  the  8.4  GHz  observations  after 
a  7  yr  interval  would  place  a  5  km  s-1  limit  on  the  velocity 
(presuming  that  no  change  in  the  scattering  disk  was  seen). 

4.  THE  SCATTERING  ENVIRONMENT 
AROUND  B 1849+005 

4.1.  Nearby  Lines  of  Sight 

By  virtue  of  the  VFA’s  large  field  of  view  at  0.33  GHz 
(«2?5),  other  highly  scattered  sources  may  be  detectable  in  the 
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Fig.  12. — Lines  of  sight  in  the  0.33  GHz  field  of  view  of  B 1849+005  that 
may  also  be  dominated  by  scattering.  The  gray  scale  shows  the  smoothed, 
continuum-subtracted  Ha  emission  from  the  Southern  H-Alpha  Sky  Survey 
Atlas  (Gaustad  et  al.  2001)  on  a  linear  scale  from  0  to  30  Raleighs.  Sources 
marked  with  crosses  appear  unresolved;  sources  marked  by  circles  appear 
resolved,  with  simple  structure  indicating  that  scattering  probably  dominates; 
and  sources  marked  by  triangles  indicate  sources  whose  structures  may  be 
more  complex  and  therefore  contain  contributions  from  both  scattering  and 
intrinsic  structure.  The  size  of  a  circle  or  triangle  is  proportional  to  the  source’s 
diameter  at  0.33  GHz.  Sources  identified  are  the  pulsars  PSR  B 1849+00, 
PSR  B 1854+00,  PSR  B 1855— 01,  the  supernova  remnants  Kes  79,  3C  391  and 
3C  392,  and  the  maser  OH  34.26+0.15.  The  source  B1849+005  is  just  above 
PSR  B 1849+00  (viz.,  Fig.  13). 

direction  of  B 1849+005.  With  only  a  single  frequency,  I 
cannot  demonstrate  unambiguously  that  scattering  dominates 
the  structure  of  other  sources  in  the  field.  However,  given  that 
B 1849+005  and  PSR  B 1849+00  (see  following  section)  are 
both  highly  scattered  and  that  PSR  B 1845— 01  also  exhibits  a 
large  amount  of  pulse  broadening  (Ramachandran  et  al.  1997), 
it  is  reasonable  to  consider  that  scattering  may  contribute  to 
the  observed  diameters  of  many  of  these  sources. 

Figure  12  shows  the  angular  diameters  for  all  sources  in  the 
0.33  GHz  field  of  view  of  B 1849+005  that  are  judged  to  be 
compact  and  likely  dominated  by  scattering.  Table  5  tabulates 
the  diameters  of  these  sources. 

If  the  structures  of  these  sources  are  dominated  by  scatter¬ 
ing,  there  is  substantial  variation  in  the  scattering  on  sub¬ 
degree  angular  scales.  Slysh  et  al.  (2001)  have  reached  a 
similar  conclusion  by  comparing  their  space  VLBI  observa¬ 
tions  of  the  maser  OH  34.26+0.15  to  scattering  measurements 
for  PSR  B 1849+00  and  B 1849+005.  They  find  the  maser  to 
have  an  angular  diameter  of  less  than  1  mas,  from  which  they 
deduce  a  scattering  measure  of  1.68  xlO-3  kpc  m-20/3,  far 
less  than  that  of  PSR  B 1849+00  or  B 1849+005,  though  in  this 
case  it  is  likely  that  the  maser  lies  in  front  of  the  strong 
scattering  region. 

Figure  12  suggests  that  the  region  of  intense,  yet  spatially 
variable  scattering  may  subtend  a  large  angle  on  the  sky.  For 
instance,  there  are  a  number  of  sources  near  (18h47?5,  +2°) 
whose  diameters  range  from  unresolved  to  1" .  Numerous 
other  cases  of  unresolved  sources  fairly  close  to  sources  with 
measurable  diameters  also  can  be  identified.  These  lines  of 
sight  may  be  similar  to  other  directions  showing  both  heavy 
and  anisotropic  scattering  and  large  changes  in  the  scattering 


strength  on  degree  size  scales  such  as  the  Cygnus  region  (Fey 
et  al.  1989;  Lazio  et  al.  1990;  Spangler  &  Cordes  1998;  Desai 
&  Fey  2001),  the  Galactic  center  (van  Langevelde  et  al.  1992; 
Frail  et  al.  1994),  and  along  the  Galactic  plane  (Dennison  et  al. 
1984;  Fey  et  al.  1991). 

4.2.  The  Line  of  Sight  to  PSR  B 1849+  00 

The  pulsar  PSR  B 1849+00  is  only  13'  from  B 1 849+005  and 
is  the  most  heavily  scattered  pulsar  known,  with  a  pulse 
broadening  time  of  0.2  seconds  at  1.4  GHz.  Lohmer  et  al. 
(2001)  show  that  between  1.4  and  2.7  GHz  the  pulse  broad¬ 
ening  of  this  pulsar  is  anomalous  in  the  sense  that  its  fre¬ 
quency  dependence  is  v~x  with  x  =  2.8^^  rather  than  the 
expected  x  «  4.  In  turn,  this  anomalous  frequency  scaling 
implies  a  density  spectral  index  of  a  =  6. 6+23  •  (For  a  review 
of  other  observations  of  this  pulsar,  see  Stanimirovic  et  al. 
2003.) 

My  observations  of  B 1849+005  are  at  frequencies  that 
bracket  those  used  by  Lohmer  et  al.  (2001)  to  observe 
PSR  B 1849+00  (Fig.  7).  Unfortunately,  the  large  uncertain¬ 
ties  on  the  angular  diameters  between  0.33  and  2.3  GHz 
preclude  definite  conclusions  about  the  possibility  of  similar 
anomalous  scattering  along  the  line  of  sight  to  B 1849+005. 
On  the  one  hand,  the  frequency  scaling  for  the  angular  di¬ 
ameter  between  1.5  and  2.3  GHz  is  v~23  and  the  frequency 
scaling  between  0.33  and  5  GHz  is  v~2 .  Both  of  these  fre¬ 
quency  scalings  (as  well  as  Fig.  7)  would  suggest  a  frequency 
dependence  of  approximately  v~2  is  a  much  better  description 
for  the  angular  broadening  than  a  dependence  of  approximately 
v~XA  (because  the  pulse  broadening  and  angular  broadening 
are  related  as  r  ~  02).  Recall,  though,  that  the  0.33  GHz  di¬ 
ameter  may  be  overestimated  because  of  contamination  from 
the  source  structure  (§  3.2),  the  1.5  GHz  diameter  is  somewhat 
smaller  than  the  resolution  of  the  array,  and  at  2.3  GHz  the 
diameter  has  been  estimated  from  only  a  single  baseline. 

Cordes  &  Lazio  (2001;  their  Figs.  2  and  3)  showed  that 
anomalous  pulse  broadening,  over  a  limited  frequency  range 
and  possibly  consistent  with  that  observed  by  Lohmer  et  al. 
(2001),  can  be  produced  by  two  scattering  screens  having 
different  scattering  strengths.  Such  a  scenario  would  produce 
anomalous  angular  broadening  as  well.  I  model  this  anoma¬ 
lous  scattering  in  a  piecewise  linear  fashion,  with  angular 
broadening  scaling  approximately  as  v~2  2  below  1 .4  GHz  and 
above  2.7  GHz,  while  between  1.4  and  2.7  GHz  it  would  scale 
as  v~XA .  This  model  produces  an  overall  frequency  scaling 
between  0.33  and  5  GHz  of  v~X  9 .  Thus,  while  Figure  7  shows 
that  the  angular  diameter  appears  to  scale  as  expected  with 
frequency,  a  modest  amount  of  anomalous  scattering  at  fre¬ 
quencies  near  1.5  GHz  cannot  be  excluded. 

A  more  stringent  constraint  on  the  anomalous  scattering  of 
PSR  B 1849+00  could  be  obtained  from  multifrequency  ob¬ 
servations  of  its  angular  diameter.  The  field  of  view  of  the  VLA 
is  sufficiently  large  that  PSR  B 1849+00  is  detected  in  both  the 
0.33  and  1.5  GHz  observations.  Unfortunately,  at  0.33  GHz  it 
is  unresolved,  while  at  1.5  GHz  it  is  sufficiently  near  the  edge 
of  the  primary  beam  that  it  is  difficult  to  obtain  a  reliable 
angular-diameter  estimate.  VLBI  observations  are  under  way 
(and  will  be  reported  elsewhere)  to  observe  B 1849+005  and 
PSR  B 1 849+00  simultaneously  near  1 .5  GHz  in  an  effort  to  see 
if  the  pulsar’s  angular  diameter  appears  anomalous. 

Nonetheless,  by  combining  the  various  observations  for 
B 1 849+005  and  PSR  B 1 849+00,  it  may  be  possible  to  constrain 
the  properties  of  the  regions  responsible  for  the  excess  scat¬ 
tering.  The  scattering  toward  B 1849+005  and  PSR  B 1849+00 
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TABLE  5 

Other  Possibly  Scattered  Sources  Near  B 1849+005 


Right  Ascension 

Declination 

Galactic  Longitude 
(deg) 

Galactic  Latitude 
(deg) 

Diameter 

(arcsec) 

(5) 

Structure 

(6) 

18  35  10.617 . 

01  45  52.86 

32.650 

4.422 

? 

18  37  18.258 . 

-01  50  26.68 

29.680 

2.302 

3.8 

18  41  44.789 . 

-01  52  52.35 

30.152 

1.297 

5.1 

18  43  58.583 . 

01  53  27.68 

33.766 

2.522 

3.7 

18  44  07.675 . 

00  33  17.10 

32.592 

1.879 

0.6 

18  44  50.864 . 

01  37  15.89 

33.624 

2.205 

3.4 

18  45  17.789 . 

00  07  42.94 

32.117 

1.308 

5.2 

? 

18  46  03.708 . 

00  03  39.24 

32.265 

1.168 

3.3 

18  46  05.681 . 

00  25  38.53 

32.703 

1.384 

3.1 

18  46  14.538 . 

-01  19  40.83 

31.158 

0.550 

5.0 

18  46  23.171 . 

00  18  28.74 

32.630 

1.264 

18  46  31.286 . 

00  09  08.35 

32.507 

1.163 

18  46  39.510 . 

00  16  23.36 

32.630 

1.188 

18  46  41.266 . 

00  55  14.98 

31.571 

0.637 

18  46  53.543 . 

00  41  02.80 

33.023 

1.323 

2.6 

? 

18  47  09.999 . 

01  41  55.15 

33.957 

1.724 

18  47  43.843 . 

01  54  33.63 

34.209 

1.695 

3.0 

18  47  55.076 . 

01  22  19.10 

33.752 

1.409 

7.0 

18  48  23.488 . 

-01  23  58.55 

31.339 

0.039 

4.2 

18  48  25.281 . 

00  46  36.09 

31.897 

0.317 

5.7 

18  48  34.187 . 

00  36  18.25 

32.067 

0.362 

5.5 

18  48  45.681 . 

01  52  55.68 

34.302 

1.453 

3.1 

18  49  12.392 . 

02  11  32.05 

34.629 

1.496 

4.1 

18  49  33.461 . 

01  54  29.11 

34.416 

1.288 

2.6 

18  50  10.398 . 

00  20  04.02 

33.086 

0.434 

2.7 

18  50  20.140 . 

00  49  13.91 

33.537 

0.619 

2.6 

18  50  40.088 . 

03  57  33.61 

36.370 

1.975 

2.4 

18  51  12.962 . 

00  27  37.89 

33.317 

0.259 

3.5 

18  51  52.819 . 

00  40  01.42 

32.389 

-0.403 

3.9 

18  52  12.978 . 

02  22  52.19 

35.140 

0.912 

5.1 

? 

18  52  27.364 . 

00  32  00.70 

33.523 

0.017 

18  53  23.434 . 

00  21  37.42 

32.834 

-0.599 

1.5 

18  53  36.039 . 

00  27  14.06 

32.775 

-0.688 

1.9 

18  53  37.212 . 

02  24  56.68 

35.331 

0.616 

6.0 

? 

18  53  43.689 . 

00  19  56.05 

33.489 

-0.358 

18  54  34.997 . 

02  27  51.16 

35.484 

0.424 

18  54  59.951 . 

02  15  39.94 

35.350 

0.239 

6.4 

18  56  09.016 . 

00  02  04.43 

33.500 

-1.033 

7.5 

18  56  27.229 . 

01  35  57.48 

34.927 

-0.386 

3.6 

18  56  31.711  . 

00  17  27.47 

33.771 

-1.000 

3.2 

18  56  43.059 . 

00  32  50.93 

34.021 

-0.925 

18  57  16.674 . 

00  11  21.39 

33.766 

-1.213 

2.3 

18  57  21.043 . 

00  45  28.03 

34.281 

-0.970 

18  57  49.365 . 

00  51  17.91 

34.421 

-1.031 

1.8 

18  58  07.936 . 

01  25  49.86 

34.969 

-0.837 

1.3 

18  58  44.966 . 

00  55  59.52 

34.596 

-1.201 

6.6 

? 

18  59  08.028 . 

00  41  48.99 

34.430 

-1.394 

4.0 

Notes. — Units  of  right  ascension  are  hours,  minutes,  and  seconds,  and  units  of  declination  are  degrees,  arcminutes,  and 
arcseconds  in  J2000  coordinates.  Col.  (5):  unresolved  sources  noted  with  ellipses  (. . .).  Col.  (6):  question  mark  indicates  whether 
a  visual  inspection  of  the  source  image  suggests  a  contribution  from  intrinsic  structure  to  the  measured  angular  diameter. 


is  a  factor  of  5  or  more  in  excess  of  what  is  predicted  by  the 
large-scale  components  (a  thick  disk,  thin  disk,  and  spiral  arms) 
of  the  Cordes-Lazio  model  (and  the  model  by  Taylor  &  Cordes 
1993).  In  order  to  reproduce  the  observed  scattering  toward 
PSR  B1 849+00,  Cordes  &  Lazio  (2003)  introduced  a  “clump” 
or  “cloud”  of  electrons  along  the  line  of  sight  to  it.  They  require 
the  clump  to  have  an  internal  electron  density  ne  =  10  cm-3  and 
a  fluctuation  parameter  F  =  200.  The  fluctuation  parameter  F 
(see  also  Cordes  et  al.  1991)  is  a  combination  of  the  fractional 
variation  of  the  electron  density  internal  to  the  clump  and  the 


outer  scale.  Larger  values  of  F  mean  that  a  clump  is  more 
efficient  at  scattering.  (Lor  a  clump,  the  fluctuation  parameter 
and  the  scattering  measure  are  proportional,  SM  oc  Fn2eAr, 
where  A r  is  the  size  of  the  clump.)  They  placed  the  clump  at 
a  distance  of  7.4  kpc  and  assumed  its  size  to  be  A r  =  20  pc. 

The  power  of  combining  the  measurements  is  twofold. 
Lirst,  the  magnitude  of  angular  broadening  is  strongly  de¬ 
pendent  upon  the  distribution  of  scattering  material  along  the 
line  of  sight,  with  material  close  to  the  observer  contributing 
more  strongly  (van  Langevelde  et  al.  1992;  Lazio  &  Cordes 
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1998).  Consider  a  Galactic  source  at  a  distance  D  with  angular 
broadening  #Gai  and  an  extragalactic  source  along  a  similar  line 
of  sight  with  angular  broadening  0xgai-  If  the  scattering  along 
the  line  of  sight  is  dominated  by  a  single  screen  at  a  (radial) 
distance  A  from  the  Galactic  source,  then  (van  Langevelde 
etal.  1992) 


#Gal  _  A 
$xgal  F) 


(5) 


This  relationship  implies  that,  even  though  seen  through  the 
same  scattering  screen,  a  Galactic  source  will  have  a  smaller 
scattering  diameter  than  an  extragalactic  source;  this  result 
occurs  because  the  scattering  screen  is  less  effective  at  scat¬ 
tering  the  diverging  wavefronts  from  the  Galactic  source  than 
the  planar  wavefronts  from  the  extragalactic  source. 

The  0.33  GHz  scattering  diameter  of  B 1849+005  is  4"2. 
Taking  an  upper  limit  on  the  scattering  diameter  of  PSR 
B1849+00  to  be  3"  (=1/2  the  beam  diameter),  A/D  <  0.7.  As 
PSR  B1 849+00  is  at  a  distance  of  roughly  8  kpc,  A  <  5.7  kpc, 
or  the  scattering  screen  must  be  more  than  2.3  kpc  from  the 
Sun. 

The  second  reason  for  combining  the  measurements  toward 
B 1849+005  and  PSR  B 1849+00  is  that  angular  broadening 
and  pulse  broadening  have  different  weightings  for  the  dis¬ 
tribution  of  scattering  material  along  the  line  of  sight.  If  x  is 
the  fractional  distance  between  the  pulsar  and  the  observer 
(with  x  =  0  being  at  the  pulsar  and  x  =  1  being  at  the  ob¬ 
server),  then  pulse  broadening  is  weighted  as  x(l  —  x),  while 
angular  broadening  is  weighted  as  x2. 

I  have  repeated  the  analysis  of  Cordes  &  Lazio  (2003)  but 
now  requiring  that  a  single  clump  describe  the  excess  scat¬ 
tering  toward  both  PSR  B 1849+00  and  B 1849+005.  Specifi¬ 
cally,  I  adopt  the  large-scale  components  of  the  Cordes-Lazio 
model  and  then  search  for  a  clump  of  excess  scattering  that 
can  explain  the  observed  scattering  diameter  for  B 1849+005, 
the  upper  limit  on  the  scattering  diameter  of  PSR  B 1 849+00, 
and  the  pulse  broadening  of  PSR  B 1849+00.  A  grid  search 
over  the  three  parameters  ne,  F,  and  x  was  used  to  search  for 
the  best  fit  in  a  minimum  %2  sense.  The  size  of  the  clump  is 
assumed  to  be  40  pc,  so  that  it  covers  both  lines  of  sight  even 
when  the  clump  is  near  the  pulsar. 

An  added  constraint  in  the  fitting  process  is  that  the  dis¬ 
persion  measure  of  the  pulsar  (DM  =  680  pc  cm-3)  and  the 
modeled  distance  to  the  pulsar  must  be  self-consistent.  That  is, 
increasing  the  density  within  the  clump  necessarily  reduces 
the  modeled  distance  to  the  pulsar,  which  in  turn  affects  the 
weighting  for  the  scattering  observables.  The  fitting  procedure 
was  therefore  the  following.  The  dispersion  measure  contrib¬ 
uted  by  the  clump,  <5 DM,  was  determined  for  each  set  of 
values  for  ne  and  x.  The  model  distance  for  the  pulsar  was 
determined  by  requiring  that  the  large-scale  components  of  the 
Cordes-Lazio  model  account  for  the  remainder  of  the  pulsar’s 
DM,  i.e.,  for  (DM  —  <5 DM).  This  new  distance  was  then  used 
in  determining  the  weighting  factors  for  the  modeled  scatter¬ 
ing  observables.  In  practice,  the  <5 DM  contributed  by  the 
clump  is  sufficiently  small  that  it  does  not  affect  the  modeled 
pulsar  distance  («9  kpc)  significantly. 

No  set  of  parameters  can  reproduce  all  three  observational 
constraints.  In  particular,  the  pulsar’s  pulse  broadening  is 
consistently  underestimated  (by  about  an  order  of  magnitude), 
while  the  extragalactic  source’s  angular  diameter  is  over¬ 
estimated  (by  about  a  factor  of  2-3).  As  might  be  expected, 
because  the  observational  constraints  are  all  due  to  scattering, 
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Fig.  13. — Subimage  of  Fig.  12  showing  the  lines  of  sight  to  B1849+ 
005  and  PSR  B 1849+00.  The  gray  scale  shows  the  smoothed,  continuum- 
subtracted  Ha  emission  on  a  linear  scale  from  0  to  10  Raleighs.  The  sizes  of 
the  symbols  are  proportional  to  the  scattering  diameter  (B1 849+005)  or  upper 
limit  (PSR  B 1849+00)  at  1  GHz.  For  the  scattering  diameter  of  B 1849+005,  it 
is  shown  with  the  correct  orientation  and  anisotropy. 


the  parameters  ne  and  F  are  highly  correlated.  The  best-fit 
value  for  the  distance  to  the  clump  is  4.7  kpc.  While  consistent 
with  the  crude  estimate  obtained  above  based  on  just  the 
scattering  diameters,  this  distance  for  the  clump  corresponds 
to  x^  1/2.  This  value  simply  reflects  the  value  of  x  that 
maximizes  the  modeled  pulse  broadening,  which  still  only 
poorly  reproduces  the  observed  pulse  broadening. 

If  the  excess  scattering  toward  both  B 1849+005  and 
PSR  B 1849+00  is  attributed  to  the  same  clump,  the  angular 
diameters  provide  the  most  useful  constraint:  the  clump  must 
be  at  least  2.3  kpc  distant.  Given  that  no  single  clump  repro¬ 
duces  all  of  the  observations,  it  is  quite  likely  that  different 
clumps  affect  the  two  lines  of  sight. 

Is  there  a  structure  along  the  line  of  sight  that  could  con¬ 
tribute  to  anomalous  scattering  toward  PSR  B 1849+00  but  not 
B 1849+005?  Figure  13  shows  the  Ha  emission  in  a  small 
region  around  these  two  sources  taken  from  the  Southern 
H- Alpha  Sky  Survey  Atlas  (SHASSA;  Gaustad  et  al.  2001). 
Clearly  apparent  between  the  two,  but  offset  slightly  toward 
PSR  B 1849+00,  is  a  clump  of  slightly  stronger  Ha  emission. 

This  Ha  clump  appears  to  lie  on  the  outskirts  of  the  H  n 
region  G33.418— 0.004  (Lockman  1989).  If  this  Ha  clump  is 
associated  with  the  H  n  region,  the  Ha  clump  would  represent 
only  a  small  fraction  of  the  total  Ha  luminosity  of  the  H  n 
region.  Based  on  the  recombination  line  strengths  for  the  H  n 
region,  much  of  it  must  be  heavily  absorbed  at  optical  wave¬ 
lengths.  This  line  of  sight  is  toward  the  Aquila  Rift  (Dame 
et  al.  2001),  and  the  area  to  the  south  of  the  Ha  clump  appears 
consistent  with  being  heavily  obscured  (Fig.  13). 

The  emission  measure  of  this  H  n  region  is  EM « 
7000  pc  cm-6  (Lockman  1989).  The  Ha  clump  and  the  line  of 
sight  to  PSR  B 1849+00  appear  to  be  where  the  5  GHz  con¬ 
tinuum  has  dropped  by  approximately  a  factor  of  2  (Altenhoff 
et  al.  1979).  Thus,  the  EM  contributed  by  the  H  n  region  along 
the  line  of  sight  to  the  pulsar  is  perhaps  3500  pc  cm-6.  The 
equivalent  scattering  measure  is  log  SM  =  0.81  [kpc  m_20/3], 
assuming  an  outer  scale  of  1  pc  (Cordes  et  al.  1991).  This 
level  of  scattering  is  comparable  to,  though  somewhat  below, 
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that  required  for  the  line  of  sight  to  PSR  B 1849+00.  The 
fitting  procedure  above  produces  a  clump  with  log  SM  =  1 .4, 
but  the  scattering  parameters  of  the  clump  in  the  Cordes-Lazio 
model  imply  log  SM  =  3.6. 

Although  it  seems  reasonable  to  attribute  some  fraction  of 
the  scattering  toward  PSR  B1849+00  to  G33.418— 0.004,  the 
H  ii  region  does  not  appear  to  be  capable  of  explaining  it  all. 
The  line  of  sight  to  PSR  B 1849+00  would  have  to  pierce 
several  hundred  similar  H  n  regions.  Higher-density  filaments 
or  tendrils  from  the  H  n  region  may  be  in  front  of  the  pulsar 
and  responsible  for  the  anomalous  scattering,  but  the  identi¬ 
fication  of  such  filaments  must  await  Ha  observations  of 
higher  resolution  and  sensitivity. 

Alternately,  a  larger  SM  would  be  obtained  if  the  outer  scale 
in  G33. 418— 0.004  were  l0  1(T4  pc  (/0  20  AU).  Given 

that  the  size  of  G3 3. 41 8— 0.004  is  probably  comparable  to  1  pc, 
it  is  difficult  to  understand  why  the  outer  scale  would  be  so 
much  smaller.  Moreover,  the  H  n  region  is  well  offset  from  the 
location  of  B 1849+005  (Altenhoff  et  al.  1979),  so  the  enhanced 
scattering  along  the  line  of  sight  to  the  extragalactic  source 
cannot  be  attributed  to  the  H  n  region.  As  Spangler  et  al.  (1986) 
concluded,  the  enhanced  scattering  toward  B 1849+005  (or 
PSR  B 1849+00)  cannot  be  attributed  to  any  single  object  or 
region. 

5.  CONCLUSIONS 

I  have  observed  the  extragalactic  source  B 1849+005  at 
frequencies  between  0.33  and  15  GHz  with  the  Very  Large 
Array  (VLA)  and  Very  Long  Baseline  Array  (VLB A)  and 
reanalyzed  archival  VLA  observations  between  0.33  and 
4.9  GHz.  The  source  has  a  complex  structure,  as  Ley  et  al. 
(1991)  suggested,  with  components  of  emission  both  to  the 
southwest  and  northeast  of  the  central  component.  However,  at 
least  to  1 5  GHz,  the  central  component  is  compact  and  highly 
scattered;  our  larger  frequency  range,  nearly  two  orders  of 
magnitude,  extends  the  earlier  analyses  of  Spangler  et  al. 
(1986)  and  Ley  et  al.  (1991),  who  also  found  the  source  to  be 
highly  scattered. 

I  have  analyzed  the  interferometric  visibilities  to  extract 
parameters  of  the  spectrum  of  density  fluctuations  responsible 
for  scattering  and  compared  the  scattering  of  B 1849+005  to 
other  nearby  sources.  The  density  fluctuations  are  anisotropic, 
with  an  anisotropy  of  about  1.3,  comparable  to  that  seen  for 
other  highly  scattered  sources  in  the  Galactic  disk  but  less  than 
that  toward  Galactic  center  sources.  While  the  data  are  sug¬ 
gestive  that  the  position  angle  of  the  anisotropy  changes  with 
frequency,  the  uncertainties  are  large.  Nonetheless,  it  seems 
clear  that  density  fluctuations  are  becoming  anisotropic  on 
scales  of  order  1015  cm  ( D/5  kpc).  An  inner  scale  to  the 
density  spectrum  also  appears  necessary  to  reproduce  the 
observed  visibilities.  The  evidence  for  an  inner  scale  l\  comes 
both  from  direct  fits  for  a  density  spectrum  including  an  inner 
scale  as  well  as  from  assessing  the  values  of  the  density 
spectral  index  a  for  a  power-law-only  density  spectrum.  Both 
of  these  approaches  suggest  an  inner  scale  of  roughly  a  few 
hundred  kilometers,  comparable  to  that  found  toward  other 
sources. 

The  two  observations  at  0.33  GHz  span  approximately  a 
decade.  There  is  essentially  no  change  in  the  shape  of  the 
image  over  this  interval.  A  conservative  upper  limit  on  the 
velocity  of  the  scattering  material  is  1800  km  s-1.  This  value 
is  sufficiently  large  that  it  does  not  exclude  material  from  the 
supernova  remnant  G33. 6+0.1  from  contributing  to  the  en¬ 
hanced  scattering  along  this  line  of  sight. 


In  the  0.33  GHz  image  of  the  field  around  B 1849+005, 
there  is  a  number  of  other  sources  that  might  also  be  heavily 
scattered.  As  is  seen  in  other  regions  of  intense  scattering 
(e.g.,  the  Cygnus  region  and  the  Galactic  center),  there  are 
large  changes  in  the  strength  of  scattering  on  lines  of  sight 
separated  by  a  degree  or  less.  With  only  a  single  frequency  to 
estimate  the  scattering,  not  all  of  the  large  diameters  in  this 
region  can  be  attributed  to  scattering.  Nonetheless,  many  of 
the  sources  probably  are  heavily  scattered,  and  these  other 
sources  in  the  field  could  serve  as  a  useful  starting  point  for 
future  scattering  observations  in  this  direction. 

Both  B 1849+005  and  PSR  B 1849+00,  which  are  sepa¬ 
rated  by  only  13',  show  extremely  strong  scattering,  and  PSR 
B 1 849+00  shows  an  anomalous  frequency  scaling  for  its  pulse 
broadening  between  1.4  and  2.7  GHz.  While  there  is  no  indi¬ 
cation  of  an  anomalous  frequency  scaling  for  the  angular 
broadening  of  B 1849+005  above  5  GHz,  the  observations 
between  1.4  and  2.7  GHz  (namely  those  at  1.5  and  2.3  GHz) 
are  insufficient  to  place  strong  constraints  on  the  presence  of 
any  anomalous  frequency  scaling.  Moreover,  the  observations 
were  insufficient  for  placing  useful  constraints  on  the  angular 
broadening  of  PSR  B 1849+00;  simultaneous  observations  of 
PSR  B 1849+00  and  B 1849+005  at  1.6  GHz  are  under  way 
and  will  be  reported  elsewhere. 

If  the  lines  of  sight  to  the  pulsar  and  extragalactic  source  are 
affected  by  the  same  “clump”  of  scattering  material,  it  must 
be  at  least  2.3  kpc  distant.  However,  a  detailed  attempt  to 
account  for  all  of  the  scattering  observables  toward  these 
sources  (angular  broadening  of  the  extragalactic  source,  pulse 
broadening  of  the  pulsar,  and  upper  limits  on  the  angular 
broadening  of  the  pulsar)  does  not  produce  a  self-consistent 
set  of  parameters  for  a  clump  that  can  reproduce  all  three 
measured  scattering  observables. 

There  is  a  clump  of  Ha  emission  that  lies  in  between  the  lines 
of  sight  to  PSR  B 1 849+00  and  B 1 849+005,  which  is  associated 
with  an  unobscured  part  of  the  H  n  region  G33.4 18—0.004. 
However,  unless  conditions  within  this  H  n  region  are  unusual, 
it  appears  unable  to  produce  a  sufficient  amount  of  scattering  to 
explain  the  excess  scattering  along  both  (or  even  either)  lines  of 
sight.  As  Spangler  et  al.  (1986)  concluded,  the  scattering  to¬ 
ward  either  B 1 849+005  or  PSR  B 1 849+00  cannot  be  attributed 
to  a  single  object  or  region. 
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